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Abstract 

We examine the Kosterlitz-Thouless universality class and show 
that essential scaling at this type of phase transition is not self- 
consistent unless multiplicative logarithmic corrections are included. 
In the case of specific heat these logarithmic corrections are identified 
analytically. To identify those corresponding to the susceptibility we 
set up a numerical method involving the finite-size scaling of Lee- 
Yang zeroes. We also study the density of zeroes and introduce a new 
concept called index scaling. We apply the method to the Xy-model 
and the closely related step model in two dimensions. The critical 
parameters (including logarithmic corrections) of the step model are 
compatable with those of the Xy-model indicating that both models 
belong to the same universality class. This result then raises questions 
over how a vortex binding scenario can be the driving mechanism for 
the phase transition. Furthermore, the logarithmic corrections identi- 
fied numerically by our methods of fitting are not in agreement with 
the renormalization group predictions of Kosterlitz and Thouless. 
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1 Introduction 



In lattice field theory one is interested in the phenomenon of phase transi- 
tions, where the quantity representing the length scale of relevant physics 
(correlation length or inverse mass gap) diverges. Typically, in spin mod- 
els, the transition is between phases with and without long range order. In 
models with a temperature driven phase transition these phases are distin- 
guished by negative and positive reduced temperature (t) respectively. A 
conventional phase transition is then characterised by a power-law diver- 
gence in the (infinite volume) correlation length near criticality (near t = 0); 



In two dimensional spin models with continuous symmetry group and con- 
tinuous interaction Hamiltonian, the existence of a phase with conventional 
long range order is precluded by the Mermin- Wagner theorem . Therefore 
there is no spontaneous magnetisation in 0(n)-spin models for n > 2. Phys- 
ically, the reason for this is that any long range order which would otherwise 
be present is destroyed by spin wave excitations in two dimensions. 

In such models there can however exist topological long range order. In- 
deed, in a (i-dimensional theory, if the order parameter lies in a space G, then 
topological defects of dimension p can occur if the [d — p — homotopy 
group TTd-p-i{G) is non trivial 0. Thus the two dimensional 0(2)-model 
can have point defects or vortices. 

Kosterlitz and Thouless used approximate renormalization group [RG] 
methods to show the existence of a phase transition driven by the binding of 
such vortices in the two dimensional XF-model (which is also called the 0(2) 
non-linear cr-model or the two component Heisenberg model) at finite non- 
zero temperature 0, ^ . The two dimensional Xy-model remains the generic 
model for the study of Kosterlitz-Thouless [KT] type phase transitions. 

The partition function for the classical n-component Heisenberg model 
on a lattice A of linear extent L is 



oo 



~ t 



(1.1) 



1.1 The Xy-Model 
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where 

d 

xeA 11=1 xgA 

and where (3 = 1/kT is the (reduced) inverse temperature (T being the 
temperature and k the Boltzmann constant). The n-component spins ffx 
have unit modulus and his a (reduced) external field. This is the infinite bare 
coupling limit of the 0{n) 0^-model. It is believed that the n-component 
Heisenberg model is in the same universality class as the general 0{n) 0^- 
model. 

The scenario proposed by Berezinskii and Kosterlitz and Thouless 
[|] is that at temperatures above some critical value (/3 < Pc) the vortices 
and antivortices are unbounded and serve to disorder the system. The vor- 
tex chemical potential is a relevant variable. Decreasing the temperature 
(/? Pc) causes the vortices and antivortices to bind, thereby decreasing 
their relevance as dynamical degrees of freedom. The model remains critical 
(thermodynamic functions diverge) for all (3 > (3c and the critical exponents 
are dependent on temperature. In terms of the reduced temperature, which 
is defined as 

*^l-f , (1.4) 

the leading critical behaviour of the correlation length, susceptibility and the 
specific heat is given in |^ as 

UiP) ~ e^*-^ , (1.5) 
Xoo(/3) ~ d , (1-6) 
Coo(/3) ~ + constant , (1.7) 

where for t 0^, v = 1/2, 7 = 7/// = 2 — 1]^ = 7/4 and a = a/u = —d = —2. 

The exponential scaling behaviour of ( |1.5| ) is referred to as essential scal- 
ing to distinguish it from the (more usual) power-law behaviour of ( |1.1| ). Thus 
the 'KT scenario' means a phase transition (i) driven by a vortex binding 
mechanism and (ii) exhibiting essential scaling behaviour. 

Although essential scaling is also often referred to as KT scaling, this (ex- 
ponential) scenario was known to exist in certain models before the seminal 
work of IQ. In particular, the BCSOS and F- models are special cases of the 
6-vertex model and are exactly solvable in the absence of an external field 
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[^]. This means that the even thermodynamic functions and in particular the 
free energy and the correlation length are exactly obtainable in these models. 
The latter has the form ( |1.5D (also with v = 1/2) and the singular part of 
the free energy is 

/r°^(/5)«C^ • (1.8) 

This model has not been solved in the presence of an external field and 
the odd critical exponents (7, rj and S) are not, therefore, rigorously deter- 
minable. (Neither are these critical exponents determinable from the usual 
scaling relations since they don't all hold in the case of an essential critical 
point 1,11). 

The reasoning presented by KT supporting an essential scaling scenario 
for the XF-model (see also p, 0) 0) 0) is non-rigorous and a primary concern 
of lattice (non-perturbative) field theory has been its verification. An unam- 
biguous verification of all the KT predictions has until now however proved 
elusive. Recently, however, Hasenbusch et al. [|l^ numerically matched the 
RG trajectory of the dual of the XF-model to that of the BCSOS model. 
This provides further evidence that at least the even critical exponents of the 
XF-model are those given by KT. 

Most Monte Carlo [MC] and high temperature expansion [TE] analyses of 
the XF-model have concentrated on the numerical determination of u and 
r) from (|1.5| ) and ( |1.6|) as they stand. These determinations are notoriously 
difficult not least because multi-parameter fits are required. Commonly, /3c is 
determined by fixing 1/ = 1/2 in ( |1.5| ). Typically, however, the subsequent fit 
to (|1.6| ) yields a value significantly different from the KT prediction rjc = 1/4:. 



Unconstrained four-parameter fits are hampered by slow convergence prob- 



lems due to extremely shallow valleys in the 4-parameter space and also 
fail in general to quantitatively and unambiguously confirm the KT predic- 
tions. Table |l] demonstrates this by listing some of the more recent estimates 
for and rjc from Monte Carlo [MC] and high temperature expansion [TE] 
methods. 



In a recent letter [^ we presented a very general theoretical argument 
showing that the essential scaling scenario is not self-consistent unless multi- 
plicative logarithmic corrections are included. We believe that the measured 
values of rjc in table 1 deviate from the theoretical value {rjc = 1/ 4) because 
these logarithmic corrections have been ignored in these analyses (the only 



4 



Table 1: Estimates for /?c and rjc for the Xy-model from a selection of papers. 
Those marked with an asterisk favour conventional as opposed to essential 
scaling behaviour. 



Authors 


Year 


Method 


f3c 


rjc 


Fernandez et al. [|12|] 


1986 


MC (FSS) 


1.1 


0.24(3) 


Seller et al. [ 




1988 


MC 


1.01(1) 


0.38(2) 


Gupta et al. |r 


I] 


1988 


MC 


1.114(2) 


0.34 


Wolff 


1989 


MC 


1.12(7) 


0.280(4) 


Biferale, Petronzio |T5| 


1989 


MC(FSS) 


1.112(2) 


0.243(7) 


Butera, Comi,Marchesini 


161 


1989 


TE 


1.112(6) 


0.300(25) 


Ferer, Mo |17 


1 


1990 


MC(BSRG) 


1.115(10) 


0.27(3) 


Hulsebos, Smit, Vink |18| 


1990 


MC(FSS) 


1.13 




Edwards, Goodman, Sokal |19] 


1991 


MC(FSS) 


1.130(15) 




Gupta, Bailie ^ 


1991 


MCRG(FSS) 


1.119(5) 


0.235(5) 


Butera, Comi l2l|] 


1993 


TE 


1.118(3) 


0.27(2) 


Janke, Nather 


22| 


1993 


MC(FSS) 




0.2389(6) 


Catterall, Kogut, Renken 2^] 


1993 


MC(FSS) 




0.247(1) 


Butera, Comi ^ 


m 


1994 


TE 




0.27(1) 


Schultka, Manousakis [^] 


1994 


MC(FSS) 


1.119(4) 




Hasenbusch, Marcu, Finn 


10| 


1994 


MC(BSRG) 


1.1197(5) 




Olsson 


m 


1995 


MC 


1.1209(1) 




Kim |27 


1 * 


1995 


MC 


1.108(3) 




Campostrini et al. |28| 


1995 


TE 


1.118(4) 


0.228(1) 



exception being [^). Our assertion is that it is essential to take account 



of the possibility of multiplicative logarithmic corrections in any attempt to 
numerically verify an essential scaling scenario. 

In we also presented a new method involving the finite-size scaling 
[FSS] of the first Lee- Yang zero |Q whereby the critical point f3c and the 
leading and logarithmic correction critical indices can be determined numer- 
ically from two parameter (straight line) fits. Here we wish to elaborate on 
these theoretical and numerical approaches and to extend our analysis to the 
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first ten Lee- Yang zeroes. 



1.2 The step Model 

The step model or sgn 0(2) model is closely related to the XF-model pi 



|, Its partition function is given by ( p..2D with 



^^^''^ = E E sgn(a- ■ . (1.9) 

xgA 11=1 

The question of criticality of this model has, until now, been unresolved de- 
spite several analyses based on high temperature series and on numerical 
simulation [^, The interest in the model arises from its possible mem- 
bership of the Kosterlitz-Thouless universality class. Like the Xy-model, 
the step model has a configuration space which is globally and continuously 
symmetric. Unlike the Xy-model, however, the interaction function is dis- 
continuous and the Mermin-Wagner theorem remains unproven for this 
case (see however [0). Nonetheless, it is expected that if a phase transition 
exists in the step model, it should not be to a phase with long range order 



® ©Hi. 

The KT phase transition of the XF-model is understood to be driven 
by the binding/unbinding of vortices. The energetics of vortex formation in 



the step model are, however, very different from the XF-model |3^, ^ 
Since vortices with effectively zero excitation energy can be created at all 
non-zero temperatures, the usual KT argument does not naturally lead one 
to expect such a phase transition in the step model. 

Sanchez- Velasco and Wills ||3^ presented evidence of critical behaviour 



starting at = 0.91 ±0.04. This was based on FSS of the spin susceptibil- 
ity. Since the associated critical index r7(/5^*°P) was significantly greater than 
that measured for the XY-model, it was concluded that the step and XY 
models are not in the same universality class. In this paper we present evi- 
dence that the step model is not in fact critical at that temperature. However 
it is critical at lower temperatures with a critical index Ti{(3f^^) compatible 
with the value. 

These results are consistent with the possibility that the XY- and step 
models belong to the one universality class. Since the KT vortex unbinding 
mechanism is not believed to be responsible for the phase transition of the 
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step model, the same scenario has to be questioned in the XF-model. This 
is not the first time the KT scenario has been questoned. For earher counter- 
evidence to both the physical KT picture as well as the quantitative essential 
scaling predictions see [HOl, O, E3] and |]13|, |39|. Further criticisms of the 



conventional view are found in 43 



The accuracy afforded by the Lee- Yang zeroes study and the consider- 
ation of multiplicative logarithmic corrections are essential for the present 
analysis [El . 



2 Lee— Yang Zeroes and Logarithmic Correc- 
tions to Scaling 

The subject of partition function zeroes was introduced in 1952 by Lee and 
Yang as an alternative way to understand the onset of criticality in 
statistical physics models. For a finite system the zeroes of the partition 
function are strictly complex (non-real). As L is allowed to go to infinity 
one generally expects these zeroes to condense onto smooth curves. Zeroes 
in the plane of complex external magnetic field h are refered to as Lee- Yang 
zeroes. Lee and Yang further showed that for certain Ising-type systems these 
zeroes are in fact restricted to the imaginary h axis (the Lee- Yang theorem) 
Dunlop and Newman proved that the Lee- Yang theorem holds for the 
two dimensional XF-model While no specific proof exists for the step 



model, its similarities to the XF-model lead one to anticipate a similar locus 
of zeroes. In practice, one obtains numerical confirmation of this hypothesis 
by explicit determination of the locus. 

In the symmetric phase, t > 0, the Lee- Yang zeroes lie away from the real 
/i-axis, pinching it only as t \ (in the thermodynamic limit). For systems 
obeying the Lee- Yang theorem this pinching occurs at /i = 0, prohibiting 
analytic continuation from 'Re{h) < to 'Re{h) > 0. 

2.1 Lee— Yang Zeroes 

Consider a partition function of the form 

ZL{(3,h) = ^Y.^'''''''-"' > (2-1) 
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where n is a unit vector defining the direction of the external magnetic field 
and h is a scalar parameter representing its strength. The factor A/l ensures 
that Zl{0,0) = 1. Here S and M are given by ( |1.3| ) or (|1.9|), S taking real 
values from the interval [—dL'^, dL'^] {dL'^ is the number of links in A) and 
M = ifi.M) taking the real values in [— L'^, L'^] [L'^ being the number of sites 
in A). 

Suppose now that the M-range is binned such that there are 2Nm M-bins 
of width Am = L'^/Nm- Defining the integrated spectral density p^^p{M) as 
the number of configurations having given M- values, the partition function 
for this system can be written as 

Z[''\f3,h) = '--^ pf}{kAM-M^^^)[e^^-^) , (2.2) 

where Mmax = (Nm — |)Am and where k e N and J\fi^^ is such that 
Z^j^\o,0) = 1. Thus Z^j^\(3,h) is proportional to a polynomial with real 
coefficients of order 2M^a.x/ Am in the 'fugacity' z = exp (/iAm) and as such 
has 2Mmax/AM zeroes in the complex z-plane. 

For the Ising model (for example) Am is the minimum change in the 
total magnetization upon flipping a single spin. There Am = 2. In the 
XF-model the M-range is continuous, even for a finite system. However, 
when we come to the numerical analysis of the XF-model, the Monte Carlo 
integration method can only sample a discrete subset of the infinite range of 
configurations open to the system. It thus accesses only a discrete sample of 
M-values from the range [—L'^,L'^]. If no binning is used in the numerical 
approach, then the minimum separation between two such M-values will play 
the role of Am above. In any case, the partition function is still proportional 
to a polynomial in the fugacity. The real and imaginary parts of the partition 
function are separately periodic in lm{h) with period I-k j Am- The pattern 
of Lee- Yang zeroes is itself therefore periodic in \m.{K) with the same period. 

2.2 Scaling and Corrections to scaling 

Assume that the Lee- Yang theorem holds for the binned system. In terms 
of its zeroes, 

z,(/3,L) = e'^«^^(^'^) , (2.3) 
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the partition function ( p.2| ) is 
Assume 



2Afmax/AA./ 

n {z-z,{(3,L)) 



fL{(3M = lim = lim L-'^lnZf 



(2.4) 



(2.5) 



Zeroes on the unit circle in the complex fugacity plane correspond to zeroes 



on the imaginary h axis. Following p9| , |50|, |51[], we define the density of these 
zeroes as 



(2.6) 



i=i 



Scaling in the Thermodynamic Limit Using the symmetry and peri- 
odicity of the pattern of zeroes one finds for the the free energy in the 



thermodynamic limit |46 



/oo(/?, h) = /,eg(/?, h)-2r j^^^G^iP^ e)de 



where 6'c(/5) is the Yang-Lee edge and is given by |S 

^7oo(A^) = o for -e,{(3)<e<e,{(3) 



(2.7) 



(2.8) 



The zero field susceptibility is the second derivative of the free energy with 
respect to /i, 



Xoo(/3) 



h=0 



roc 



(2.9) 



Following 1^ |5T|, |29[, the cumulative density of zeroes can be written in 
terms of the susceptibility and the Yang-Lee edge: 



MP), 



(2.10) 
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$(x) being some function of x with $(| x |< 1) = 0. 
The (singular part of the) specific heat is then 



h=0 



(2.11) 



The above formulae are quite general and hold for any model provided 
only that it obeys the Lee- Yang theorem. To proceed from here we have to 
insert the (model-specific) critical behaviour. Instead of the conventional KT 
formula ( |1.6| ) and ( |1.7| ) assume, now, the following more general behaviour 
for the zero field susceptibility, the edge and for the singular part of the 
specific heat 



Xoo{(3) 

cr^(/3) 



d(ineoo)^ 



and 



7 
a 



r = —vr 
q = —vq 



Then (|ig) and (piP give 



a-7-2A 



(In^o 



OC 




1 + 



(2.12) 
(2.13) 
(2.14) 

(2.15) 
(2.16) 



(2.17) 



Inserting the leading KT behaviour for the correlation lengthQ (|1.5| ) gives 

a = 7 + 2A , (2.18) 



and 



g = r + 2p + 2 



1 

1 + - 

V 



(2.19) 



Here we have the first indications that KT scaling behavoiur ( [1.5| ) - ( |1.7| ) 
is not self-consistent without multiplicative logarithmic corrections. If there 
are no logarithmic corrections then g = f = p = and ( |2.19| ) cannot hold. 

^Including logarithmic corrections in (1.5) so that ^ exp {bt^'^)t'^ only delivers 
extra additive corrections and has no consequence for the other critical indices. Although 
such corrections may be present in the correlation length, w remains undetermined. 
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Finite— Size Scaling The finite-size scaling (FSS) hypothesis, first formu- 
lated in 1971 by Fisher and co-workers |52[, is a relationship between the 
critical behaviour of thermodynamic quantities in the infinite volume limit 
and the size dependency of their finite volume counterparts. The general 
statement of FSS, which is expected to hold in all dimensions |j5^ is that if 
Pl{(3) is the value of some thermodynamic quantity P at inverse temperature 
(coupling) (3 measured on a system of linear extent L, then 



Poo{P) 



(2.20) 



where is the correlation length of the finite-size system. Here JFp is 

some P-dependent function of the scaling variable x = ^l{Pc) / ^oo{P) ■ Fixing 
the scaling variable, one has 

eoo(/3) ~ 2;-ia(/3c) . (2.21) 

Luck has shown that for the two dimensional XY model, C,l{i3c) is propor- 
tional to L Therefore FSS for the susceptibility and the Yang-Lee edge 



is from (gJ^) and (|2l3 



XdPc) 



L\\nLY 



(2.22) 
(2.23) 



Now, from ( p.4|) and (|2.71) the magnetic susceptibility for a finite-size system 
is 

h=0 



If the lowest lying zeroes exhibit the same FSS behaviour |55|, then 



We provide numerical justification for this assumption in Sec. ff] 
i U^ , (i]23D and {^^W^) give 



(2.25) 
Together, 



A 



P 



(2.26) 
(2.27) 
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whence 

a = -d , (2.28) 

and 

g = 2(^l + ^^ . (2.29) 

Thus the KT result for the leading specific heat index a is recovered. However 
the correction index is g = 2(z/ + 1) = 3, a non-trivial result. There are 
therefore multiplicative logarithmic corrections to the leading (KT) scaling 
behaviour of the singular part of the specific heat. 

Cr^~C'(lneoo)'~CV^ • (2.30) 

The full specific heat has another (constant) term coming from the regular 
part of the free energy and since the leading critical index is negative any 
numerical verification of ( p.30| ) is rendered very difficult if not impossible 
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Integrating ( p.30|) twice with respect to t gives the singular part of the 
free energy to be foo ~ ^oo^ as in (pT^ ). Thus we have indirectly verified 
hyperscaling for the two-dimensional XF-model and we have done this using 
the self-consistency of the essential scaling behaviour and finite-size scaling^. 

The above analytic considerations have yielded no information on the 
odd correction exponent r. The original renormalisation group analysis of 
Kosterlitz and Thouless @], in fact, implicitly contained the prediction 

r = -l/16 (2.31) 



as noted by [|2T[. Subsequent analyses have concentrated on the r = form 
of the scaling behaviour ( |2.22D and the verification that r]{(3c) = 1/4. Allton 
and Hamer [^] have conjectured that the deviation of their determination 
of 7] from 1/4 might be due to logarithmic corrections. 

With u = 1/2 and A = —2 + r]c/2 the FSS formula for the first Lee- Yang 
zero is 

e^m^L^lnLY . (2.32) 

The study of this full scaling form and the numerical determination of r is 
the subject of what follows. 



^ Hyperscaling and FSS are distinct hypotheses, the latter being the stronger. In the 
upper critical dimension {d = 4) FSS has been shown to hold where hyperscaling fails p3|. 
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3 Numerical Results 



The Monte Carlo (MC) method is a stochastic method of importance sam- 
phng and an intrinsically non-perturbative approach to the calculation of 
path integrals and expectation values. In this section we wish to report on 
an application to the Xy-model in an attempt to test independently the 
above scaling scenario. 

In recent years the development of more efficient techniques has enhanced 
the quality and practicality of certain numerical calculations in bosonic spin 
systems. These improvements include 

• the advent of non-local algorithms M, EB[ which greatly reduce critical 



slowing down, 

• the rediscovery and development of the spectral density method which 
allows extrapolation away from the simulation point ||5^, |60| , 

• the use of multihistograms to extend the coupling range over which extrap- 



olation of information is reliable 1611, |62 



The aim of the present section is to present numerical evidence for the 
existence (or otherwise) of logarithmic corrections to scaling. We shall adopt 
a self-consistent strategy. 

1. (a) Temporarily ignore the effect of logarithmic corrections and aim 

to extract the basic critical parameters Pc and rj{(3c). One verifies, 
from the scaling of Xl{P), the existence of a critical region below 
some approximately determined temperature = l/jSc and that 
the effective values of 7 near this temperature include the expected 
KT value 2-7]^ = 7/4. 

(b) Assuming this value holds at (3c-, use it to extract the so-called 
Roomany-Wyld beta-function approximant ( 3.38| ) which is then 
valid only at the critical point. 

(c) Deduce (3c from the zeroes of this function using the residual L 
(finite-size) dependence to estimate systematic errors. 

2. (a) Now use multi- histogramming techniques to study the scaling of 

Lee- Yang zeroes (with much higher precision) over a range of 
candidate (3c values around the above best estimate. 

(b) Verify that a critical region exists and that the leading behaviour 
is compatible with the KT value (A = —15/8) 
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(c) Use the quality of the scahng fits to 6^1 (/5c) of ( |2.32| ) at each can- 
didate Pc to estabhsh the possible existence of logarithmic correc- 
tions. 



3.1 Simulation Details 

A non-local algorithm based on that of Wolff fl^ was used to simulate the 
XY- and step models at zero magnetic field {h = 0) on square lattices of sizes 
L = 32, 64, 128 and 256. In the case of the step model, additional simulations 
were performed at L = 48. The values of f3 at which the simulations took 
place (Po) were evenly spaced in steps of 0.02 between 0.90 and 1.20 for the 
XF-model and with various degrees of spacing between 0.86 and 1.40 for 
the step model. At each simulation point, 100, 000 measurements of 5* and 



M were made using ( |1.3| ) and ( |1.9| ). 

In these isotropic two component models (with no external field) the 
orientation of the magnetization is not fixed in an y direction. One uses for 
the magnetization its magnitude M = + where and My are the 
1- and 2-components of the magnetization in the two dimensional internal 
spin space . 



3.2 Approximate Critical /3 Determination 

The expected finite-size scaling behaviour of critical systems can be used to 
locate critical points and determine critical parameters of an infinite system. 
A convenient method for accomplishing this |63, p5[ is to extract an approx- 



imation to the Callan-Symanzik beta function from the size dependance 
(in lattice units) of some physical quantity. 

Consider the measurement of a physical quantity Q on a lattice of linear 
extent C = La where a is the lattice spacing in physical units. The mea- 
surement of Q on the lattice at bare (dimensionless) coupling (3 will result 
in some number expressed in lattice units: Ql{i3)- Renormalisation may be 



viewed |Q as varying the scale unit a while maintaining physical quantities 
such as C and Q fixed at their physical (bulk) values. Thus, the lattice gran- 
ularity L = C/a varies when a does. Fixed physics then requires that the 
bare coupling /? is correspondingly tuned and the measure of its response is 
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just the Callan Symanzik [CS] beta function 



B{P) = -f^p-'f-, /i = - . (3.33) 
a/i a 

Roomany-Wyld approximant Roomany and Wyld ||65| showed how, for 
Q = A4 = 1/^ (the mass-gap), one can readily extract an estimate of B{P) 
from finite lattice measurements of Ql{P) = fniiP) = 1/^l(/9)- The gener- 
alisation to other physical quantities was made by Hamer and Irving p^. 

To find a preliminary estimate for /3c (step 1(a) above), we ignore logarith- 
mic corrections and write FSS for the magnetic susceptibilty (from (|2.22|) ) 
as 

^approx _ ^2 _ ^3 3^^ 

We consider the physical, dimensionless combination 

V = xC"- = Lha^ = L^lW) . (3.35) 
The corresponding CS equation is 

dV dV 

dV = -^dp + -^dL = (3.36) 
op oL 

where L is treated as a continuous variable, for the present. The correspond- 
ing beta function is then deduced as 

^ .d\nP ,d\nP 

Numerical approximations to the derivatives are then applied as circum- 
stances dictate. Using lattice sizes L and L' to estimate the logarithmic 
derivative one obtains the generalised RW approximant 

R , <PQ-HQLit)/Q,{t)]/\n{L'/L) 

BQ{t,L,L} = 1 a 1 r/^ r+\r^ mi ("^"^^^ 

iftHQL'{t)QL{t)] 

Classical numerical interpolation may be used to estimate the derivative ap- 
pearing in the denominator. 
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Application to the XY— Model The present application is to the deter- 
mination of the critical point using finite lattice data for the susceptibility 
alone (ignoring multiplicative logarithmic corrections). Following through 
the above arguments for the expected beta function critical behaviour, one 
finds 

(3.39) 

This is independent of 7 but the construction of the numerical approximant 
(|3.38|) requires its knowledge (0q = 7/1^). There are two further complica- 
tions 

• the presence of logarithmic corrections has been ignored up to this 
point 

• there are several arguments to show that 7 = 2 — is weakly dependent 
on t. The KT prediction is that, at t = {P = Pc) , V = \- 

Thus the RW beta function will only be even approximately valid at the 
location of the critical point. Since it should vanish there, this is good enough 
to locate the critical point subject to the assumed value of the index 7. 



Results The above has been carried out for the XY-model using lattice 
sizes L = 32, 64, 128 and 256. Specifically, the effective 7(T, L) has been 
extracted from 3 pairs of consecutive L values using the scaling implied by 
(|3.34|) . The difference between effective values was monitored and seen to 
vanish above (3 = 1/T ^ 1.1. The average values were also montitored and 
seen to pass through 1.75 at the same place. Statistical errors were estimated 
assuming that those of the input XiiP) data are Gaussian and independent. 

Next, treating 7 = 1.75 as constant in the neighborhood of (3 = 1.1, the 
RW approximant to the beta function was calculated for consecutive pairs 
of available L values. As expected, these vanished at stable values of (3 as 
shown in Fig. |I|. The estimates of critical (3, which we now denote (3bw, were 
determined by numerical interpolation and statistical errors determined as 
above. The values are shown in Table ^ 



The spread of values with the choice of L, L' used to determine them is 
clearly very small and can be used to estimate a systematic error (within 
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Table 2: Preliminary /3rw estimate from RW approximant zeroes. 





XY Model 


step Model 


L/L' 




stat. err. 


^step 


stat. err. 


32/48 






1.1977 


0.0095 


48/64 






1.2089 


0.0094 


32/64 


1.1040 


0.0014 






64/128 


1.1045 


0.0014 


1.2104 


0.0052 


128/256 


1.1064 


0.0011 


1.2184 


0.0068 



the approximation of ignoring multiplicative logarithmic corrections). We 
estimate 

/5rw = 1.106 ± .002 ± .002. (3.40) 

where the errors are respectively statistical and systematic. 

We emphasise that we ignore logarithmic corrections in this section and 
/3rw is therefore only a preliminary or approximate estimate for (5c. 

Application to the Step Model A similar analysis was performed for the 
step model using lattice sizes L = 32,48,64, 128 and 256. The RW approxi- 
mant results are also shown in Table 2 and the corresponding (3^^ estimates 
in Fig. lb. Although the data are somewhat noisier, the behaviour of the 
step model is strikingly similar to that of the XY-model. The prehminary 
estimate of critical P is 

= 1.21 ± 0.01 ± 0.01 . (3.41) 

3.3 Lee— Yang Zeroes and Numerical Test of Scaling 
Scenario 
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3.3.1 Determination of Lee— Yang Zeroes 

When the external field is complex [h = h^. + ihi), the partition function 
( |1.2|) can be rewritten 

Zl{P, K + ihi) = ReZLiP, K + ihi) + ilmZL{P, K + ihi) , (3.42) 

where 

ReZL{(3, K + ih^) = ^ E e^^+^^^ cos {h.M) , 

= Z LiP,K) {cos {h,M))(,,hr, (3.43) 

and 

lmZL{(3,hr + th,) = ^^e^^+'^'-*^sin(/i,M), 

{^x} 

= ZLiP,K){smih,M))f,,h,.. (3.44) 

According to the Lee- Yang theorem the zeroes are on the imaginary /i-axis 
{hr = 0) where (sin (/ijM))^ vanishes and so the Lee- Yang zeroes are 
simply the zeroes of 

ReZLiP, ihi) oc (cos (/iiM))^,o ■ (3.45) 

Thus the Lee- Yang zeroes are easily found and at no stage is a simulation 
with complex action involved. 

Errors were calculated using the straightforward bootstrap method where 
the data for each Po are resampled A^boot = 100, times (with replacement) 
leading to A^boot estimates for 6i, from which the variance and bias can be 
calculated. In fact, to take account of any correlation present in the 100,000 
values of S and M at each /Sq, the following method of error determination 
was (also) used. 

At each /3o the 100,000 data are split into Asput subsets each containing 
100, 000/Aspiit data. The bootstrap method is then applied to each of these 
A'spiit subsets using A'boot = 50 in each case. If the data is strongly correlated 
or if the error is strongly dependent on whatever correlation is present, the 
dependence of the error on A'sput should be visible. In fact, plots of the relative 
error in the position of the zero against A^sput at each f5o and at each L reveal 
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little or no Ngpht dependency. Thus the data are highly uncorrelated and the 
error estimates reliable. Moreover, within the range of /3o-values studied, 
there appears to be only a weak dependence of the relative errors on /3o (the 
relative errors decrease slightly with increasing (3 as the zero approaches the 
real axis) and none discernable on L. For the XF-model these relative errors 
are 0.0003 for the lowest lying zeroes, increasing to 0.0020 for the 15th zeroes 
while for the step model they increase from 0.0004 to 0.0030. 



Multihistograms and the Spectral Density Method In order to be 
able to find the Lee- Yang zeroes at temperatures other than the simulation 
points /3o, we employ a multihistogram reweighting technique (the spectral 



density method |6T], ^) which accomodates extrapolation between 

values. In view of the large amounts of data now involved it is necessary 
to introduce binning. For the spectral density p^^\S,M) Sec. 2.1 we bin the 
raw histograms in a 500 x 100 array. The results for the Lee- Yang zeroes 
turn out to be very stable with respect to the bin sizes. Altering the number 
and size of bins has only a tiny effect on the position of the zero (well within 
the eventual errors). 



3.3.2 The XF-Model 

Scaling of Lowest Lee— Yang Zeroes We begin our analysis with the first 
Lee- Yang zeroes at the preliminary critical coupling estimate /3rw — 1.11. 
In Fig. ^ we plot the logarithm of the position of the first Lee- Yang zero 
against the logarithm of the lattice size L, at /3 = 1.11. In the absence 
of any corrections, the slope should give the leading power-law exponent 
A. In fact at f3 = 1.11 the slope is Xes = —1.8778(2). The deviation from 
the KT value of —15/8 = —1.875 can be attributed to (a) the presence of 
logarithmic corrections and/or (b) the possibility that (3 = 1.11 is not in 
fact the critical point. As a preliminary investigation of (a) we identify the 
correction exponent r in ( |2.32| ), by plotting ln(^iL'^^/^) against In(lnL) in 
Fig. ||. A straight line is identified. Its slope is —0.012(1) giving evidence 
for a non-zero value of r, albeit not in agreement with the renormalization 
group predictions of —1/16 = —0.0625 from 0. 

The above estimate /5rw — l-H for the critical coupling came from an 
analysis in which the possibility of the existence of logarithmic corrections 
was ignored. At this point, we must allow for the possibility of the true 
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critical coupling jSc being different than /3rw- To systematically investigate 
the extent to which the above conclusions depend on the measurement of (3^ 
we now study the above FSS (fits for Aefr) and corrections to FSS (fits for r) 
as a function of the assumed critical beta. 

In Fig. 4a we plot Aefr coming from fits to the leading FSS behaviour 
as a function of 13. Fig. 4b gives the corresponding per degree of free- 
dom [x^/dof] for these fits. Not all fits yield an acceptable x^/dof. To be 
conservative, we take 'acceptable' to mean 

xVdof < 2 , (3.46) 

which corresponds to a minimum confidence level of 14%. For Fig. ^ we have 
acceptable values of x^/dof only for [3 ^ 1.107. The figure indicates that 
FSS sets in above this value. I.e., the system remains critical for (3 ^ 1.107. 

This range of critical [3 values corresponds to Aeg ^ —1.8761(2). This 
range of critical (3 values does not include that for which the KT value of 
— 15/8 = —1.875 is recovered (Fig. ^ gives Aeg ^ —1.875 at /3 ~ 1.105 with 
X^/dof ~ 3.4). Thus we can conclude that ignoring logarithmic corrections 
gives deviation from the KT value of —15/8 = —1.875 for A. 

Next we investigate how our conclusions regarding the correction expo- 
nent r depend on the critical coupling. Fig. ^ gives the results from fits to r 
from (|2.32| ) and the corresponding x^/dof. we see that acceptable values are 
possible only for 1.107 ^ [3c ^ 1.119. The r = solution again corresponds 
to (3c — 1.105 and x^/dof ~ 3.4. Similary, we find that a fit with r = —1/16 
would correspond to [3c — 1.138 and x^/dof ~ 13. 

Therefore we have numerical evidence challenging the detailed quantita- 
tive predictions of Kosteritz and Thouless. In summary, the results of the 
analysis of the first Lee- Yang zeroes of the XF-model is that assuming the 
KT value = 1/4, we find non-zero logarithmic corrections to scaling and a 
corresponding estimate of the critical temperature: 

r = -0.02(1), /3c = 1.113(6). (3.47) 



KT Versus Power Law Scaling Scenario The KT essential scaling sce- 
nario was called into question in |T3[ where numerical evidence in support 
of conventional power law scaling behaviour in the 2'(10) clock model was 
given. While this model is expected to be in the same universality class at 
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the Xy-model (same critical indices) it may have a different critical coupling 
f3c hindering a comparison with our numerical results. However a recent pa- 
per by Kim also claims to favour a convential as opposed to KT phase 



transition in the XF-model. 

It is extremely difficult to distinguish between KT and power law scaling 



on the basis of numerical methods alone |T9[. Power law scaling for the 
thermodynamic functions (with additive power law corrections) would lead 
to (|2.28|) being replaced by Josephson's law 

a = 2-vd , (3.48) 

and the FSS of the zeroes ( p.32| ) by 

^i(/3e) =aiL^(l + a2^"') , (3.49) 

for some Oi, 02 and 03. If z/ > 1 the singular part of the specific heat decreases 
with lattice size as in the KT case ( p.30|) . Thus Josephson's law may not 



be a numerically useful criterion for distinguishing between the two scaling 
scenarios. 

The leading FSS of 61 in (|3.49|) is the same as for the KT case ( p.32| ) 



and ignoring corrections leads to the same conclusions regarding Aefr- How- 
ever, as seen above, it is essential to include corrections to correctly identify 
Pc- Ignoring corrections would lead us to conclusions similar to The 
numerical analysis of ||2^ applies to and x- There, rjl/S) is measured by 



fitting to the FSS formula xl ~ ^^^^ (ignoring multiplicative logarithmic 
corrections). Good fits (x^/dof < 0.5) were reported for (3 > 1.105. Using 

X = -2 + ri/2 , (3.50) 

the corresponding X^s can be found. These are listed and compared to our 
own results (when corrections are ignored) in Table |^. The agreement is 
impressive. 

In the critical point is located by finding a coupling, /?p, where xl ~ 
l2~r,W) j^Qi^g ^ ^ ^/4^ 'j^Yie resuh is 1.105 < Pp < 1.111 with 'most 
probably' /5p = 1.1062(6) (T^ = 0.9040(5)). This compares well to our result 
r = (no corrections) at /5 ~ 1.105. Using this as the critical coupling, 



reports that fits to ^ in (|1.5| ) give an unacceptable x /dof. In the light of 
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Table 3: Comparison of our results for X^s (ignoring multiplicative logarith- 
mic corrections) to [27]. 



T 


(3=1/T 


1] from 27| 


Xcs from [27] 


Our Acff 


Our x^/dof 


0.910 


1.099 


0.2600(10) 


-1.8700(10) 


-1.8712(2) 


9.0 


0.905 


1.105 


0.2520(05) 


-1.8740(05) 


-1.8750(2) 


3.4 


0.900 


1.111 


0.2450(05) 


-1.8770(05) 


-1.8778(2) 


0.8 


0.890 


1.1236 


0.2310(10) 


-1.8840(10) 


-1.8838(2) 


0.2 



our analysis this is hardly surprising. The system is still in fact in the high 
temperature phase and criticality has not been reached - jSp is not the true 
critical coupling. 



Examination of power law scaling in [^ in the form of a four parameter 
fit to ^oo = Af^i^l + bt^) yields a reasonable x^/dof. We attempted a three 
parameter fit to ( p.49|) and found acceptable x^/dof over a range of values 
for (3 but believe this to be attributable to the fact that it is easy to make 
a three parameter fit to only 4 data points. So while we cannot logically 
exclude a power law scaling scenario on the basis of our numerical results 
neither can the (modified) KT scaling scenario be excluded on the basis of 



The lesson here is that it is essential to include consideration of logarith- 
mic corrections to KT scaling in order not to be led to the wrong value of jSc 
upon which conclusions are highly sensitive. 

3.3.3 The Step Model 

Scaling of Lowest Lee— Yang Zeroes The analysis began with a rough 
search for the leading critical behaviour predicted by ( |2.32| ) . An independent 
test was also made using the (less accurate) susceptibility data and (|2.22|) . 
Both methods indicated critical behaviour setting in for (3 > 1.2. 

In Fig. 6 we display the results for the effective exponent A^g' as a 
function of (3 together with the corresponding x^/dof. Acceptable fits are 
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only found for ^ 1.185. We note that the corresponding values of A^g' 
( < - 1.8715(2)) include that (-15/8 = -1.875) corresponding to the KT 
prediction. 

Again, the evidence for critical behaviour is the existence of a range of 
acceptable chi-squared values for a linear fit. In view of the similarity of A^g^ 
to the expected KT value, we now proceed to test the further hypothesis 
that the step model is in fact in the same universality class as the XY- 
model. We assume the behaviour ( |2.32| ) with A^*'^'' = —15/8 and search for 
corrections in Fig. ^ We find acceptable fits for 1.195 < jS^^'^^ < 1.295 with 
0.009 > r > -0.034. 

The range of acceptable r values includes that found earlier for the XY- 
model and that corresponding to no logarithmic corrections (r = 0). Again, 
this range excludes the prediction r = —1/16 coming from the approximate 
renormalisation group treatment of the XY-mode\ 0. Thus we conclude 
that the present data are compatible with the step model being in the same 
universality class as the XF-model. We do not, however, exclude other 
possibilities. 

In summary, 

XF-Model : 1.107 < pc ^ 1-119 , -0.01 > r > - 0.03 , (3.51) 
Step Model : 1.195 < /3f^P < 1.295 , 0.01 > r''''^ > - 0.03 . (3.52) 

3.3.4 Higher Order Corrections 

We have shown that a FSS analysis ignoring multiplicative logarithmic cor- 
rections leads to deviation from the KT prediction —15/8 for A or, equiva- 
lently, 1/4 for rj{(3c)- When multiplicative logarithmic corrections are taken 
into account the results are compatible with ?7(/?c) = 1/4 but, now, the lead- 
ing logarithm exponent r is not that given by KT. 

In the same spirit, it is conceivable that numerical measurement of these 
leading logarithmic corrections may become compatible with the KT RG 
predictions (r = —1/16) if higher order corrections are accounted for. 

Keeping higher order (additive) corrections in the RG treatment of the 
magnetic susceptibility modifies the right hand side of (|2.12| ) by H, ^ 
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The FSS behaviour of the Yang-Lee edge is now expected to be (from ( |2.25|) ) 



OM) = a,L' (In LY 1^1 + a2^-^^ . (3.54) 

Unbiased four parameter fits are hampered by the convergence problems 
mentioned in section 1.1. However, if we again accept the leading scaling 
behaviour A = —15/8, two parameter fits to ai and a2 in ( |3.54| ) yield un- 
acceptable x^/dof if r is —1/16 in both the XY- and step models. In fact, 
X^/dof ranges from 13 to 2 as /3 decreases from 1.14 to 1.10 in the Xy-model 
and from 9 to 2.5 as f3 decreases from 1.40 to 1.17 in the step model, indi- 
cating that such fits may become acceptable for very low f3 in both cases. 
(Note that in the case of the XF-model such low P is outside the range of 
/^c^measurements in table 1). On the other hand, if we fix r = —0.02 we find 
X^/dof ^ 2 over the whole /5-range in both models. 



4 Analysis of Higher Index Zeroes 



In the calculation leading to ( p.25| ) we assumed that the low index zeroes 



have the same FSS behaviour. We now wish to (a) check that assumption, 
(b) redo the whole analysis for each index j = 2, . . . , 10 separately and (c) 
introduce a new concept we call 'index scaling' and apply it to these first ten 
zeroes. In fact we can determine up to fifteen zeroes for each lattice size. 
However the errors for the highest index zeroes become too large to allow 
reliable analysis. 



4.1 Analyses of Individual Higher Index Zeroes 



We repeat the analysis of (|2.32|) for each index separately for j = 2, . . . , 10 



and for both the XY- and step models. The results of the analyses for 
Aeff are sumarized in Table ^. For the XF-model the analysis applied to 
individual higher index zeroes (j ^ 4) appears to accomodate the RG result 
A = —1.875. However higher index zeroes are less accurately determined than 
lower index ones and the j = 1 results remain the most accurate. 

Accepting A = —1.875 and looking for corrections from applications of 
the r-analysis to each index j separately gives the results summarized in 
Table ^ Acceptable fits are obtained only in the (3 range given and lead to 
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Table 4: The results of the Aeg-analysis for the j-index zeroes for the XY- 
and step models: The criterion ( p.46 ) holds for each j for (3 above and A 
below the values listed. 





XF-Model 


ste 


p Model 


j 


P 


AcfT 




Aeff 


01 


1.107 


-1. 


8761(2) 


1.185 


-1. 


8715(2) 


02 


1.107 


-1. 


8761(2) 


1.185 


-1. 


8716(2) 


03 


1.106 


-1. 


8754(2) 


1.180 


-1. 


8716(2) 


04 


1.106 


-1. 


8753(2) 


1.185 


-1. 


8717(2) 


05 


1.104 


-1. 


8741(2) 


1.190 


-1. 


8725(2) 


06 


1.103 


-1. 


8735(2) 


1.200 


-1. 


8739(2) 


07 


1.103 


-1. 


8735(2) 


1.215 


-1. 


8757(2) 


08 


1.103 


-1. 


8735(2) 


1.215 


-1. 


8757(2) 


09 


1.100 


-1. 


8717(2) 


1.195 


-1. 


8735(2) 


10 


1.100 


-1. 


8717(2) 


1.185 


-1. 


8722(2) 



the listed ranges for r. These (3 and r ranges are graphically represented in 
Fig. ^ (for the Xy-model) and Fig. 9 (for the step model). Although in both 
models the higher index zeroes individually allow for r = 0, at no stage (up to 
j = 10) do we find results compatable with the RG result —1/16 = —0.0625. 
So the higher index zeroes analysed individually yield the same information 
with decreasing accuracy as j increases providing further strong evidence of 
the validity of our earlier conclusions. 

Finally, analyses of higher order (additive) corrections for the individual 
higher index zeroes yield the same qualitative results as those presented in 
sec. 3.3.4. 

4.2 Analysis of the Density of Zeroes (Index Scaling) 

The last section demonstrates the success of FSS applied to each zero sep- 
arately. We can now be confident that the lowest lying zeroes scale in the 
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Table 5: The Results of the r- Analysis for the j-Index Zeroes for the XY- 
and step models: The criterion (|3.46|) holds for each j for the following j3 
ranges and corresponding r range. 





Xy-Model 


step Model 


3 




r 




r 


01 


1.107- 


-1.119 


-0.005(1) 


- -0.030(1) 


1.195- 


-1.295 


0.009(1) 


- -0.034(1) 


02 


1.106- 


-1.120 


-0.002(1) 


- -0.032(1) 


1.195- 


-1.295 


0.009(1) 


- -0.034(1) 


03 


1.106- 


-1.120 


-0.002(1) 


- -0.031(1) 


1.195- 


-1.295 


0.008(1) 


- -0.034(1) 


04 


1.105- 


-1.119 


0.001(1) 


- -0.029(1) 


1.195- 


-1.295 


0.008(1) 


- -0.034(1) 


05 


1.104- 


-1.117 


0.004(1) 


--0.024(1) 


1.200- 


-1.295 


0.005(1) 


- -0.034(1) 


06 


1.104- 


-1.116 


0.004(1) 


--0.022(1) 


1.205- 


-1.295 


0.002(1) 


- -0.034(1) 


07 


1.104- 


-1.116 


0.004(1) 


--0.022(1) 


1.210- 


-1.300 


-0.001(1) 


- -0.035(1) 


08 


1.103- 


-1.116 


0.007(1) 


--0.021(1) 


1.210- 


-1.300 


0.001(1) 


- -0.035(1) 


09 


1.102- 


-1.119 


0.009(1) 


--0.026(1) 


1.200- 


-1.310 


0.004(1) 


- -0.038(1) 


10 


1.100- 


-1.122 


0.014(1) 


- -0.031(1) 


1.190- 


-1.350 


0.009(1) 


- -0.047(1) 



same way with lattice size L and that this scaling behaviour is 

e,{(3c)-^L\lnLy . (4.55) 

In this paragraph we wish to report an attempt to combine the scaling be- 
haviour of the various indexed zeroes. 

The density of zeroes for the finite (binned) system is given in ( |2.(j| ) as a 
sum over delta functions, 

N 

gL{(3,e)=L-'J2^i(^-(^0^L)) , (4.56) 

where N is the total number of zeroes. The corresponding cumulative density 
function G is a step function with 

GL{f3,e) = {- if ^G(^„%+i) . (4.57) 
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This function can be smeared by replacing the delta functions of ( [4.56| ) by 
Gaussians 



1 ^ 1 



1 f e-0jii3,L) 

'2\ <Tj(L) 



(4.58) 



where <Jj{L) is some appropriate quantity giving the spread of the Gaussian 
around the j^^ zero. The cumulative density function is then a sum of error 
functions 



If (Jj{L) is small enough then 



2j-l 

2N 



In the thermodynamic limit ( p.lOp gives 



(4.59) 



(4.60) 



MP), 



The corresponding FSS formula is 



(4.61) 



Gl{(3c.0)^L-^^ 



Comparing with ( ^.60| ) and using N oc L'^, 



<|.-i($(l)(2j-l)) 



The density of zeroes in the thermodynamic limit is 



dd 



(4.62) 



(4.63) 



(4.64) 



In a numerical high temperature and high field study, Kortman and Griffiths 
showed that, for Ising models in the high temperature phase, this density 
exhibits power law behaviour as a function of the distance of 9 away from 
the edge g§. 



1 



e 

W) 



(4.65) 
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At /3 = /5c, 0" = 1/5 where S is the usual magnetization index. This index 
governs the response of the magnetization to the presence of a weak external 
magnetic field, 

M ~ /i^ . (4.66) 

One therefore is led to expect 

$(x)~(l-x)^ , (4.67) 

or, 

<!>-\y)-l^y^ . (4.68) 

There are additive higher order corrections to formulae ( [4.61[ )-( p:.66| ). Now, 
([4.63| ) gives 



^ = ci + C2(2j-1)4t , (4.69) 

for some Ci and C2 which incorporate the higher order corrections. Using the 
standard scaling relation 

^ = ^^^ , (4.70) 
we can express the leading behaviour in terms of A 

^ = Ci + C2(2j-l)-^ . (4.71) 

Again, one expects there to be additive volume dependent corrections to this 
formula. 

Nonetheless, to leading order the L-dependence drops out of the ratio of 
the j^^ to I''* zeroes. This opens up the possibility of a new type of scaling 
technique - the use of the dependency of this ratio on the index j to find 
critical indices. If successful, this method could be applied in parallel to 
traditional FSS. The advantage would be that one doesn't need enormous 
lattices to eliminate finite-size effects, nor does one need a large range of 
lattice sizes to exploit FSS. In cases where only one lattice size is available, 
index scaling may, in principle, still yield a result. 



We can check the form (^4.69|) in one dimension. From the exact solution 
to the one dimensional Ising model with periodic boundary conditions, the 
Lee- Yang zeroes are given by |]30[ 



cos {29j) = -u^ + (1 - u^) cos (^ ^-^^ ^ TT^ , (4.72) 
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where u = exp (—2/3). In one dimension the critical point is at zero temper- 
ature or Uc = 0. Then the LY zeroes are at 



e,m=±\^7r . (4.73) 

Then, 

| = 2j-l , (4.74) 



which is of the form ( [4. 691) with 6 = oo, the correct value for 6 in one 



dimension. There is no exact solution to the Ising model in the presence of 
an external field in two or more dimensions. 



The XF-Model In Fig. 10 we test (^4.71|) for the two dimensional XY- 



model at (3c = 1.115. The plot involves ten LY zeroes for each of the four 
lattice sizes - forty data points in all. As expected, no L-dependence is 
visible. Fitting to 

^ = Ci + C2(2j-l)^« (4.75) 

gives, however, C3 = 0.9827(3), 5% away from the expected value —X/d = 
0.9375 and a rather poor x^/dof (8.7). The line in Fig. 10 is this best fit to 



In an early paper on the scaling of partition function zeroes ||5^ it was 
suggested that the scaling variable is in fact j/L'^. When we plot hj{L) 
against j/L'^ we find the data does not collapse onto a smooth curve. A later 
paper showed that this dependency only holds for high index zeroes [p^]. 
Indeed our result 

'2j 



= ( ^ ) (4-76) 



asymptotically approaches that of at high j. However, the poor x^/dof 



for the fit to the form ( 4.75| ) indicates that while it might be nearer the truth 



than m, { ^77^ ) is also not the full storyQ . 

To see how C3 depends on the j-indices used, we make fits to (|4.75|) for 
j = 2 — 5, 3 — 6 and so on to 7—10. We have a three parameter fit to 4x4 = 16 
data points in each case. The fits are now of good quality (x^/dof < 1) and 
the results for C3 are summarised in Table 13. The j-dependence indicates that 



Other arguments for this functional form are contained in 
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Table 6: The Results of the index scaling Analysis small ranges of j. 



J 


C3 


2-5 


0.991(1) 


3-6 


0.986(1) 


4-7 


0.981(1) 


5-8 


0.977(1) 


6-9 


0.973(1) 


7-10 


0.969(2) 



([4.75| ) can only be true asymptotically and encourages further investigation 
into the the full scaling form. A similar picture is found for the step model. 



5 Conclusions 



We have presented a very general (model independent) theoretical method 
to test the self-consistency of the scaling behaviour of odd and even thermo- 
dynamic functions. Application of this method to the KT scaling predictions 
for the two dimensional XF-model reveals that multiplicative logarithmic 
corrections cannot be ignored. In the case of specific heat these logarith- 
mic corrections were identified analytically. The corrections corresponding 
to the magnetic susceptibility were identified numerically from the scaling 
behaviour of the lowest Lee- Yang zeroes. 

The theoretical and numerical techniques we have developed are fully in- 
dependent of the physics of the phase transition. The conventional picture 
is that the XF-model phase transition is driven by a vortex binding mech- 
anism This is in contrast to the customary picture of the step model 

Here, while the configuration space is the same as 
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1 35, m, 38 



that of the XF-model, vortex formation is not believed to be energetically 
favourable given the discontinuous nature of the interaction function. If the 
XF-model phase transition is truely driven by a vortex interaction mecha- 
nism and if a phase transition exists at all in the step model then, barring a 
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remarkable coincidence, they are expected to belong to different universality 
classes (different scaling behaviour). 

Application of our numerical techniques to the two dimensional step 
model reveals the critical parameters (including logarithmic corrections) to 
be compatable with those of the XF-model. This is a strong indicaton that 
both models do in fact belong to the same universality class. This conclusion 
is reinforced by analysis of higher index zeroes. 

This result raises questions on how the vortex binding scenario can be 
the driving mechanism for the phase transition in the XY-model. Further- 
more, the multiplicative logarithmic corrections identified numerically by the 
analysis of up to the first ten Lee- Yang zeroes are not in agreement with the 
renormalization group predictions of Kosterlitz and Thouless. 



Note Since the preprint version of this paper first appeared, other authors 
have reanalysed their data to look for logarithmic corrections in the two 



dimensional XF-model. In particular, Janke |]70[ has applied FSS to the 
Villain version and found r = —0.027(1) in agreement with our result. High 
temperature measurements |]70| give r = +0.0560(17), close to Patrascioiu 
and Seller's result r = +0.077(46) l7| and r = +0.042(5) to +0.05(2) from 
Campostrini et al. ||7^ . 



We wish to thank W. Janke, E. Klepfish and CB. Lang for stimulating 
discussions. 
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Figure 1: The Roomany-Wyld beta function approximants for (a) the XY- 
model and (b) the step model. The symbols □, o, x and • 
are in order of increasing lattice size. 
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Figure 3: Corrections to FSS of the First Lee-Yang Zeroes for the Xy-model 
at 1.11. 
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Figure 4: Test of the Expected Leading Critical Behaviour for the XY 
model: (a) The effective exponent Aeg (slope of a straight hne fit) and (b) 
the corresponding x^/dof versus (3. 
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Figure 5: Logarithmic corrections for the Xy-modeh (a) the logarithmic 
correction exponent r to the Yang-Lee edge is shown as a function of the 
assumed critical coupling /3c and (b) the corresponding x^/dof. 
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Figure 6: Test of the Expected Leading Critical Behaviour for the step modeh 
(a) The effective exponent Xfg'^ (slope of a straight hne fit) and (b) the 
corresponding x^/dof versus /3. 
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Figure 7: Logarithmic corrections for the step modeh (a) the logarithmic 
correction exponent r^*''^ to the Yang-Lee edge is shown as a function of the 
assumed critical couphng and (b) the corresponding x^/dof. 
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Figure 8: Minimum and maximum acceptable (a) j3 and (b) correction ex- 
ponent r from the r-analysis with criterion (|3.46|) for the XF-model. 
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Figure 9: Minimum and maximum acceptable (a) j3 and (b) correction ex- 
ponent r from the r-analysis with criterion (|3.46|) for the step modeL 
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Figure 10: Index scaling in the XF-model at (3c = 1.115. The errors are 
smaller than the symbols. Each symbol is in fact an overlap of four data 
points coming from the four lattice sizes demonstrating excellent L-coUapse. 



46 



